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Abstract. In this paper we discuss various possibilities of using X-ray obser- 
vations to gain information about the large-scale structure of the Universe. After 
reviewing briefly the current status of these investigations we explore different 
ways of making progress in this field, using deep surveys, large area surveys and 
X-ray background observations. 

1 Introduction 

X-ray emission in the Universe arises in intense gravity environments. At high 
galactic latitudes, active galactic nuclei (AGN) and other emission line galaxies 
dominate the source counts at all explored fluxes, with galaxy clusters being the 
second most abundant source class. 

Recent RO SAT deep surveys ([Q], [^]) have shown that most of the soft X- 
ray volume emissivity in the Universe arises at redshifts z > 1 — 2 (Q). The AGN 
and star formation rate per unit volume follow a remarkably similar evolution 
rate in the Universe (|^) and therefore they can both be used as tracers of the 
evolution of large-scale structure in the Universe. Using AGN and star-forming 
galaxies as tracers of cosmic inhomogeneities is most sensitive to intermediate 
redshifts (z ~ 2), providing a critical link between cosmic microwave background 
studies (which map the z ~ 1000 Universe) and local galaxy surveys ~ 0). 

In this paper we briefiy review the current status of the use of X-ray obser- 
vations towards the study of large-scale structure. More details are presented in 
ps]]. Then we explore possibilities of making qualitative progress in this field by 
carrying out different types of X-ray surveys. 

2 What we know so far 

At galactic latitudes | 6 | > 20° the contribution from the Galaxy to the X-ray sky 
is small: less than 10% of the X-ray background above 2 keV is due to galactic 
emission, absorption is negligible above this photon energy and a census of X-ray 
sources down to any fiux limit exhibits less than 10-20% of galactic stars. Ob- 
servations of the X-ray background at high galactic latitudes and photon energies 
above 2 keV can therefore be used to map the extragalactic X-ray sky. 



2.1 The isotropy of the X-ray background 



The all-sky distribution of the X-ray background for cosmological purposes has 
been best mapped by the HEAO-1 mission. A galactic anisotropy dominates the 
large-scale anisotropy, but this can be modelled out ([^). A dipole contribution is 
detected in the X-ray sky, in rough alignment with the direction of our motion with 
respect to the Cosmic Microwave Background frame (Q, ^). The amplitude of 
this dipole accounts for both the kinematical effect of our motion (the Compton- 
Getting effect) and the excess X-ray emissivity associated with the structures 
which are pulling us. These two effects are expected to be of the same order (^) 
and the analysis done in Q shows this to be the case. However, in an analysis of 
the ROSAT all-sky data at lower photon energies (which have the disadvantage 
of a larger contamination from the Galaxy) Plionis & Georgantopoulos find a 
dipole several times larger than the expected kinematical dipole. The difference 
between both results might be partly affected by the elimination of X-ray bright 
clusters in the Scharf et al analysis, as clusters are known to be a largely biased 
population ([0). The bias parameter derived from the XRB dipole is large {bx ~ 
3-6). 

Treyer et al have analyzed higher order multipoles of the HEAO-1 A2 
X-ray background. The discrete nature of the XRB contributes a constant term 
to all multipoles which scales as oc S'cu^, where Scut is the minimum flux at which 
sources have been excised. Treyer et al detect a signal growing towards lower-order 
multipoles which is consistent with a gravitational collapse picture, as predicted 
by [^. Excluding the dipole, this analysis yields a moderate bias parameter for 
the X-ray sources {bx ~ 1 — 2). 

On smaller (a few degrees) angular scales, probing linear scales of hundreds 
of Mpc, the 'excess fluctuations' technique has been used often in the analysis of 
the XRB. The way this works is by modelling the distribution of XRB intensities 
on a given angular scale in terms of confusion noise, plus a contribution coming 
from source clustering (^]). These studies have yielded so far only upper limits 
for the excess fluctuations: < 2% on scales of 5° x 5° (0) and < 4% on scales 
1° X 2° (H). We discuss later what is the expected signal and how it could be 
measured. 

Yet on smaller (a few arcmin) angular scales, which probe the galaxy-galaxy 
clustering scale, data from X-ray imaging telescopes has been used. The autocor- 
relation function of the XRB on these scales should refiect the clustering of high 
redshift X-ray sources in the nonlinear regime (jl^). Soltan et al jl^ have found 
a strong positive detection for angular separations 0.3-20° which is, however, dif- 
ficult to interpret as both the Galaxy and the Local Supercluster could contribute 
to this. 



2.2 Clustering of X-ray selected AGN 



Studying the clustering of X-ray selected AGN is likely to be the most direct way 
to map the structure of the X-ray sky. At soft X-ray energies this requires fairly 
deep surveys (going below ^ 10~^* ergcm"'^ s~^) as otherwise very few objects at 
z > 1, where most of the X-ray volume emissivity is produced, would be sampled. 

Carrera et al have analysed a set of 'pencil beam' medium and deep 
ROSAT images containing 200 X-ray selected AGN, sampling a redshift inter- 
val 2 ~ — 2. The net result is the detection of X-ray selected AGN cluster- 
ing which is relatively weak (the 3D correlation length is ro < 5/i~^Mpc, for 
h = i/o/(100 kms~^ Mpc~^)) and strongly evolving with redshift (faster than 
comoving). At much brighter flux limits Ml] used the ROSAT All Sky Survey 
Sources to derive a 2D correlation function that, when translated to 3D with an 
appropriate catalogue depth, yields a larger correlation length (ro ~ 6/i~^ Mpc). 

2.3 Do X-rays trace mass? 

In we compile various measurements of the bias parameter for X-ray sources 
and in particular for X-ray selected AGN and the XRB. The bias parameter is 
likely to be redshift dependent. For a simple model where all objects form at the 
same early redshift, Fry finds bx{z) — bx{0) + z(bx{0) ~ 1), which implies 
that at high z the bias parameter could be large. 

The other efi'ect that comes into play, especially when using the XRB, is that 
at low redshift clusters become more numerous and their imprint in the local 
XRB features becomes more important. As clusters are a strongly biased source 
population (|^| estimate bx{0) ~ 4) it is not surprising that the amplitude of the 
XRB dipole calls for large bias factors, but higher order multipoles (sensitive to 
more distant sources) do not. 

Within present knowledge and uncertainties, the bias factor for the AGN as 
the dominant X-ray source population, appears to take moderate values bx = 1 — 2 
at low to intermediate redshifts. Indeed at higher redshifts the AGN population 
might be more strongly biased. 



3 Deep hard X-ray surveys 

Deep surveys, particularly at hard photon energies, are a key ingredient to forth- 
coming studies af the large-scale structure of the X-ray Universe. Currently pop- 
ular models for the X-ray background assume a population of AGN with a dis- 
tribution of absorbing columns (|2l], [^, [^), where most of the X-ray en ergy 
produced by accreting black holes is absorbed and re-radiated in the infrared ( |24| , 
Eq]). Several claims have been made that the absorbed AGN population evolves 
differently than the unabsorbed one (feql, pOjl, pi|). As most of the energy con- 



tent in the XRB resides at 30 keV, it is crucial to explore harder photon energies 
than previously achieved with the ROSAT deep surveys. 




Figure 1: Simulation of a 2-10 keV deep XMM EPIC pn image of 350 ks 

XMM is the most sensitive X-ray observatory to survey the X-ray sky at 
photon energies above 2 keV. Although its point-spread-function is significantly 
worst than that of Chandra, at energies above 2 keV both instruments are photon- 
starved and then the much larger collecting area of XMM will make it more 
efficient. We have carried out extensive simulations of XMM EPIC observations 
at various depths and found that the deepest planned XMM observations (PVCal, 
GTO and AO-1) reaching 350-400 ks will not be confusion noise limited. Figure 
1 shows the resulting image in the 2-10 keV band of a simulation of a 350 ks 
XMM EPIC-pn exposure in a blank field (using the standard model |2^]). As 
the EPIC field of view is 30' in diameter, we expect to find ~ 300 sources in 
such a pointing once vignetting has been corrected for. Most of these sources are 
expected to lie at redshifts 2 > 1 — 2, from which the X-ray volume emissivity in 
hard X-rays will be derived and compared with the one, assumed so far, obtained 
with ROSAT for soft X-ray photons. 

4 X-ray background surveys 

The measurement of large-scale structure in the X-ray Universe does not nec- 
essarily require individually resolving all sources in very large areas of the sky 
down to very faint fluxes. If the X-ray volume emissivity as a function of z can 
be derived from deep surveys, fluctuation analyses of the XRB can also be used 



Figure 2: Simulation of a 4-12 keV X-ray background map (ignoring the 
galaxy), as observed with a 1 collimator with a field of view of Ideg^, 
scanning the sky for 6 months 100% efficiency. 

(p7[|). If the X-ray volume emissivity peaks at some intermediate redshift as it 
does in soft X-rays (z ~ 1 — 2), then for a fixed angular scale the XRB fluctuations 
are related almost uniquely to the value of the power spectrum of the inhomo- 
geneities in the Universe at a single comoving wavenumber. The scales to be 
probed by XMM (from a few to a few tens of arcmin) will provide a measurement 
of the k ~ (0.1 — l)ftMpc~^ regime. The power spectrum is expected to peak at 
k ~ 0.05/iMpc^^, which corresponds to an angular scale of ~ Ideg. 

In we argue that to detect the excess fluctuations of the XRB due to 
source clustering for a beam size of Ideg^, a large fraction of the sky needs to 
be surveyed. To prove that this is feasible, we have carried out simulations of 
hard X-ray source populations over the whole sky with a simple clustering model 
for the sources and measured XRB intensities (details in [Q). These intensities 
are then 'measured' with a proportional counter of Im'^ effective area during one 
complete 6-month scan of the sky at 100% efhciency and including stable particle 
background in a manner similar to the Ginga LAC observations. 

Figure 2 shows one of these simulations where the clustering has been mod- 
eled with a gaussian correlation function with comoving evolution. Ignoring data 
within I 6 |< 20°, Figure 3 shows the histograms for the XRB intensities in 3 cases: 
absence of clustering, linear clustering evolution and comoving clustering evolu- 
tion {bx = 1 in all cases). The distributions are clearly distinguishable, and the 
excess fluctuations can be determined with a very high accuracy. Indeed with a 
significantly smaller collecting area and similar circumstances, excess fluctuations 
can still be detected, but measured with larger statistical uncertainties. 

One such survey will also benefit other approaches to measure large-scale struc- 
ture in the Universe, particularly the multipole expansion, especially if a sensitive 
source survey could also be carried out. In this way, those 1 deg^ regions where 
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Figure 3: Intensity distributions for the map in Fig 2. (ignoring | 6 |< 20°) 
for no clustering (continuous line), linear evolution of clustering (dashed 
line) and comoving clustering evolution (dotted line). 

sources above a given flux are detected could be masked out for the multipole 
analysis, as discussed in [ p^ . 

5 Large-area surveys 

Direct measurements of the large-scale structure of the Universe at redshifts z ~ 
1 — 2 via X-ray observations require surveying areas of hundreds of square degrees 
to a sufficient depth. Using galaxy clusters as tracers of large-scale structure of 
the Universe presents the difficulty of the faintness of most of these objects beyond 
redshifts z ~ 1. Even XMM will require a large amount of time to do a sensible 
mapping of galaxy clusters out to these redshifts. 

Using AGN has the advantage that they have strong positive evolution up to 
z ^ \ — 2. In order to reach the redshifts where most of the X-ray emissivity is 
produced z ~ 1 — 2, AGN surveys have to go down to, at least, a 2-10 keV flux 
^ lO"^"* erg cm^^ s^^ . Figure 4 illustrates the redshift distribution for different 
flux limits assuming the fesl model. The advantage of X-ray AGN surveys over 
similar optical work (e.g. the Sloan Digital Sky Survey) is that with hard X- 
rays the absorbed AGNs can also be used up to earlier times if they evolve more 
strongly than the unabsorbed broad-line objects. 

Mapping a lOOdeg^ contiguous area of the sky with the XMM EPIC cameras, 
which have the largest field of view among all operating X-ray facilities, to a depth 
of ~ 20fes (needed to get reliable detections at > 10~^* ergcm"'^ s~^) will require 
10 Msec of effective exposure time with XMM. This will collect ~ 10^ sources. 




Figure 4: Redshift distribution of AGN for different 2-10 keV flux limits. 



A more efficient way to carry out that project is by means of a dedicated 
mission with a wide-field X-ray telescope. The Panoram-X mission proposed in 
ji^ ], would cover the whole sky to a depth ~ 10~^* erg cm~^ s~^. Of course, 
finding redshifts for a reasonable fraction of the several tens of millions of sources 
to be discovered by such mission is simply impossible. A full analysis of the 
spatial distribution of X-ray selected AGN would then require the modelling of 
the 2D distribution from these maps, using redshift distributions from the deep 
hard X-ray surveys. 



6 Outlook 

X-ray cosmology is just in its infancy. Basic questions such as what is the bias 
parameter of different classes of X-ray sources are still partly unanswered. How- 
ever the fact that most of the X-ray sky is dominated by extragalactic sources of 
which AGN are the major component make the X-ray sky especially suited for 
cosmology at intermediate redshifts. 

Making quantitative progress in this field requires not only proper use of exist- 
ing or planned observatory-type facilities (i.e., Chandra and XMM), but probably 
also dedicated missions to survey all (or most of) the sky. X-ray cosmology is now 
in a position to make specific predictions for the structure of the X-ray universe 
(once hard X-ray surveys have been carried out with Chandra and XMM). These 
surveys can then be designed and optimized to obtain detections and precise mea- 
surements of the large-scale structure of the Universe at redshifts z ^ 1 — 2. That 



would really be a major boost for cosmology at intermediate redshifts. 
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